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The non-Markovian collision integral is obtained on the base of the Kadanoff- 
Baym equations for Green functions in a form with allowance for small retar- 
dation effects. The collisional relaxation times and damping width of the giant 
isovector dipole resonances in nuclear matter are calculated. For an infinite 
Fermi liquid the dependence of the relaxation times on the collective vibration 
frequency and the temperature corresponds to the Landau's prescription. 

1 Introduction 

The dissipative properties of many-body systems, specifically transport coefficients for vis- 
cosity as well as the damping of the collective excitations, are strongly governed by the 
interparticle collisions. In semi-classical approach these collisions are described by collision 
integral in kinetic equation. It allows us to involve a viscosity into motion equations similar 
to the phenomenological Navier-Stockes ones |l]]-f|. An advantage of such approach is a 
conceptual clarity and a possibility to use many results from a general macroscopic physics. 
For correct description of the collision relaxation rates in systems with a fast variation of 



the mean field, it is necessary to incorporate the memory (retardation) effects |0]-|T^]. We 



investigate the form of the non-Markovian linearized collision integral with retardation using 
the Kadanoff-Baym equations [O, [TH| for Green functions. This method enables immediately 



to take into account self-consistent mean field. 

The non-Markovian collision term of the semiclassical linearized Vlasov-Landau equa- 
tion was obtained with the use of the Kadanoff-Baym equations for the Green functions 



in Refs. [11], [LJ]. The kinetic equation was considered in one-component Fermi liquid with 
the use of the linear approximation on deviation from overall equilibrium with a nonequi- 
librium distribution function Sf = f(p,r,t) — /(e) oc exp(— iut), where f(p,r,t) is an exact 
distribution function and / is an equilibrium distribution function at equilibrium single- 
particle energy e. An expression for linearized collision integral was obtained in the Born 
approximation on two-body collisions in the following general form: 

J(p,f,t) = J {1) (p,f,t) + J {2) (p,r,t), (1) 
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where the first component corresponds to variation of the distribution function and second 
one is connected with variation 5U = e(p,f,t) — e of mean field with e(p,f,t) for actual 
single-particle energy. 

One can see that this relationship does not agree with general form of the linearized 
Landau- Vlasov equation in Fermi liquid in the limiting case of the states with week time- 
dependence of nonequilibrium component of distribution function, when dSf/dt — > 0. In- 
deed, the linearized Landau- Vlasov equation in quasi-homogeneous systems can be presented 
in the following form: 

9Sf p f t d - 

— + —{p-—)8f = J, (2) 
at m or 

where 5f is a linear deviation of the exact distribution function f(p, r, t) from the distribution 
function of the equilibrium shape /(e) evaluated with actual single-particle energy e = e+SU, 

Sf = Sf-^SU=Sf-f{e). (3) 

The derivative d8f/dt can be omitted in Eq. (0) for slightly time-dependent states and 
the collision integral should be a functional, of the Sf, 

J = *(6j)=*.0.IL)6f, (4) 
m or 

in order to the Landau- Vlasov equation would be fulfilled in this case. The relationships ([!]) 
and @ are generally in contradictions and derivation of the form of non-Markovian collision 
integral from Refs. |TT], ITJ] can be revised. 



It was pointed in Refs. || that collisional integrals of the type J = J(5f), 



are general form of the collision integrals between quasiparticles in Fermi liquid without 
retardation and they lead to a local equilibrium state described by distribution function /;. e , 
= /(e = e + 5U). Due to this we use in the following this terminology. 

In this contribution an expression for non-Markovian collision integral of the linearized 
Vlasov-Landau transport equation from Refs. [II], [l4j] is modified in Sect. 2 for the case of 



slightly time-dependent states, i.e. in a form which allows for reaching the local equilibrium 
in system. 

In Sect. 3 calculations of the relaxation times and damping width of the collective vibra- 
tion in nuclear matter with the use collision integral with retardation effect are presented. 



2 Non-Markovian linearized collision integral within 
semiclassical approach 

In order to obtain linearized Vlasov-Landau equation with collision integral we use the mixed 
{p, f} Weyl-Wigner representation for single-particle Green functions (correlation functions) 
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G>(1,1'),G<(1,1'): 

-,i ->' 

dr ' exp(--pr ')g < (t,t;r + y.r- —); ^ 
f' = f[ — r v ,r = |(fi - fi/), 

where f(p,r,t) is one-body distribution function and Green functions are determined byfll7. 



m 



(6) 



g>(i,i , ) = -^(^(i)^ + (i , )) ) «i>*;, 

G<(1, 10 =<<* + (l0 *(!)>, ti<A- 

Here, and ^(1) are the operators of creation and annihilation of a fermion; the 

symbol 1 includes both space and time variables, namely, 1 = {fi,t} (we omit isotopic and 
spin variables); and the expectation values in Eq.(|6]) are calculated for the ground state or 
for the ensemble of initial states if a temperature of system is not zero. 

The Green functions G > (1, 1 ) and G < (1, 1 ) satisfy the Kadanoff-Baym equation 



G><(1, 1 ) = /><(!, 1 ) - 7><(1, 1 ), 
G><(l,l') = Jt<(l,l')-J}<{l,l'), 

where the single-particle potentials U(l), U(l') are determined by the relations 

C/(1)G><(1,1') = J df 2 E (f 1 ,f 2 ,t 1 )G' >< (f 2) t 1 ;l / ), 
f/(l')G >< (l,l') = J df 2 G' >< (l;f 2 ,t 1 0S Q (f 2 ,f 1 ^ 1 

with 

So(n,r 2 ,t) = -i5(n - r 2 ) J drv(\n - r\)G < (r,t 1 ;r,ti) 

+iv (|fi - r 2 \)G < (f 1 ,t 1 ;f 2 ,t 1 ), 

and v(\fi — f 2 |) for two body potential. 

The functions Jf^ and are correlation integrals of the form 

1 /■*! 
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^(1,1') 

J 2 ><(1,1') 



1 /■* 



f'dh [£>(l,2)-£<(l,2)]G><(2,l'). 

^ 2 S >< (1,2) [G > (2,l')-G' < (2,l')" 
Pdt 2 [G > (1,2)-G < (1,2)]S><(2,1') 

d^G^M) [S > (2,l')-S < (2,l') 
3 



n Jt 
1 /■*! 
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These correlation integrals have the general form with allowance for retardation effect. It is 
assumed as usually that interaction between particles start at the time to = —oo. 

In order to obtain the Landau- Vlasov equation the following suggestions are adopted 
[|TT| , [14|| : A) Born approximation for two-body scattering; B) the time variation of the 
nonequilibrium distribution function 5f oc exp(— iut) is taken periodic one with real fre- 
quency during all interval of the time changing (—00 < t' < t); C) linear approximation 
on deviation of one-body Green functions from their equilibrium values are used; D) Fermi 
system is considered as quasi-homogeneous one in coordinate space. 

The linearized Vlasov-Landau equation has the form (Q) and can be presented as [11], [14 



j t f{p,r,t) + {ej} = J{p,r,t). (11) 

Here, 

d d d d 
^ ' ^ dp df^ Of dp^ 

are the Poisson brackets and 

2 

e(p,r,t) = ^ + U{p,r,t) (12) 

is the classical energy of a particle in the mean field U(p, r, t) = e + 5U(p, r, t), where 
expression for SU can be expressed in terms of the Landau interaction amplitude F(p,p ) 
for two-body collision scattering matrix: 

g r dp' 



w = ti(W F(s?)w ' r1i) ' (13) 

where Np = 2pp m*/(g ir 2 ft 3 ), pp is the Fermi momentum, m* is the effective mass of 
nucleon and g is the spin degeneracy factor. 

The linearized collision integral has form ([!]) with (see, Eqs.(42),(43) and (45), (46) of 
Ref.0 for details) 

J«\p,r,t) = 2 / ^^jf W(fi,a,A,fl)g(Aj?) B®(p,r,t). (14) 

Here, W({pi}) = (da / dQ)4(2irh) 3 / m 2 is the probability of two-body collisions with the initial 
momenta p\ = p,f>2 and final ones P3,Pa] da/dfl is in-medium differential cross-section (in 
Born approximation); 

B W (p, f, t) = J2 Sf h (t) dQ ( { p ) [5 + (Ae + hu)+64Ae-h cu)) , 

k=i °h (15) 

E( 2 )(p,r,t)=g({/J) A( ^ (t)) {[^(Ae+^)-^ + (A^]-[(54A6-M-^4A^)]}, 
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where f k = f{pk, t)\ 9Q({fj})/df k are the derivatives of the is the Pauli blocking factor Q 
with respect to f k , 

QUIA) = (i - A)(i - h)hh - hW - h)(i - />)• 

The Q = e(pi,f) and 8Uj are the equilibrium single-particle energy and the variation of the 
mean field for particle with momentum pi respectively; Ae = Z\ + e 2 — e 3 — e 4 , A(5U) = 
5Ui + 5U 2 — 8U 3 — SU4, Ap = f>i + p 2 — P3 — P4. The equilibrium distribution function 
f k = f(p k ,r) depends on the equilibrium single-particle energy e k = e(p k ,r): fk = /(efe)- It 
equals the Fermi function evaluated at the temperature T, f(e k ) = 1/[1 + exp((<ffc — ji)/T)}. 

Note that the generalized functions 5 + , <5_ appearing in Eq. fllS]) include also integral 
contribution, 



5+(x) = — [° dr e~ iXT = — = -8(x)-—V(-), 6(x) = 6*Jx), (17) 



where S(x) is the delta function and the symbol V denotes the principal value of integral 
contribution. The integral terms of the 5±, corresponding to virtual transitions, are usually 
rejected in the J because they assumed to be included by renormalizing the interactions 
between particles ||. This corresponds to substitution of the S(x)/2 for 5± in Eq. (|15|), i.e., 
to taking into account real transitions with conservation of energy. We will consider below 
only these transitions. 

Now we will modify the expression for quantity B^K We present the nonequilibrium 
component 5f of the distribution function in the form 

5f(p 3 ,r,t) = -u(p 3 ,f,t)^-. (18) 

Then the quantity can be written as 

B^\p,f,t) = -\ £ u k dQ ({f j}) [5(Ae + hco) + 5(Ae-hu;)} 

= Iai/Q({£})^ [8(Ae + hu) + 8(Ae-hu)} - SB®, 
where Au = z/ x + z/ 2 — ^3 — 1/4, v k = v{p k , f, t) and 

SBW = l^^i^QiifMK^+n^ + siAe-nu)}} 

2 k=1 de k 



(20) 



The first component in the Eq.(l20|) determines a probability flux of colliding particles 
which is connected with possibility of variation of the energy e k when the values of other 



+ \jl Q({f j W(^ + hu>) + 8{Ae-hu)]^. 

2 k=l 0e k 
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energies (e^ and hu) are fixed. This term should be equal zero because of fixing the total 
energy in initial or final states and therefore it does not contribute to the total number of 
the collisions Af: _ 

poo 

jV(p) = deJ(p,e), p = p/p. (21) 
Jo 

A relative dynamical component of the distribution function is slowly dependent on 
energy and it can be considered (at least for low temperatures T < e F ) as a function of 
the momentum direction rather than of the momentum: = v(pk,r,t) = z/(p^, ep, f, t). 
Therefore the second component in the Eq.(^) is also negligible and the term SB^ should 
be omitted from the Eq.(fLj), 5B^ = 0. The expressions for B^ take the form 

B« (p, f, t) = ±Av Qdf^-JL [5(Ae + hco) + 8{Ae -hu)], 

2 dhu (22) 

B^(p,r,t) = lg({/ J }) A( ^ {t)) {[5(Ae+hu)-8(Ae)]-[8(Ae-hu)-6(Ae)]}- 

The shift in energy Ae by hu in the arguments of the 5-functions of the expressions (p]), 
(0) ! (H) f° r the collision integral agrees with the interpretation of the collisions in the pres- 



ence of the collective excitations proposed by Landau pOfl . According to this interpretation, 
high-frequency oscillations in Fermi liquid can be considered as phonons, that are absorbed 
and created at the two-particle collisions. 

As it was discussed in the introduction, the form of expressions [l]), (0), fl22|) for collision 
integral is not correct in general case and an additional modification of these expressions 
are needed. The incorrectness results from approximations which were made at the kinetic 
equation derivation. 

Because of assumptions both on undamped behaviour of the distribution function 
throughout all time interval (— oo < t' < t — > oo) and on small magnitude of two-body 
interaction, one can expect overestimation of retardation effects in the foregoing expression 
for collision integral. It means that from physical point of view this relationship should be 
fulfilled only in the case of small u. 

Therefore, we replace the derivatives of the form d5(Ae + Tioj)/dTiuj and d 5{Ae — 
huo)/dhuo in of the Eq.([22p by the finite differences (S(Ae + hu) - 8{Ae))/%u and 
(5(Ae) — 5(Ae — hu))/hu, respectively. Then we combine the resulting expression together 
with contribution B^ arising from mean-field variation and finally obtain the linearized 
collision integral for Fermi liquid in the following form 

J(p,r,t) = J W{{pi])5{Ap)AxQ j- . (23) 

Here, Ax = Xi + X2 — X3 ~ X^'i Xk = x(Pk, t, t) is a function determining relative deviation 
of distribution function from local equilibrium state, 5f, (|3|): 

5f = Sf- ^=5U = f(p, r, t) - f Le . = ~X% X = v + 5U. (24) 
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With the use of the algebraic relation 

[(1 - - h)hh - A/aCl - /s)(l - U) exp (^j]5(Ae ± hu) = 0, (25) 



the Eq. ( P3|) can be presented as 
J (p, f , t) = 

I d ^)^ Wm})6{A ^ Axfj2 ^ " /s)(1 " h) ^ huJ ^ ~ $ (-^' T )] > ^ 

where $(hw,T) = 5{Ae + hu))[exp(-hu/T) - l}/2hu. 

The collision integral of the form ( p3| ) or ( |26| ) depends on the variation 5f, J = J(5f). 
It was mentioned in the introduction this behaviour is in line with general properties of 
the Vlasov-Landau equation in the Fermi-liquid 0, |3| at ddf/dt = 0. It provides driving 
distribution function towards its local equilibrium value. This behaviour is in line with 
general properties of the Vlasov-Landau equation in the Fermi-liquid @, [| at dSf /dt = 0. 
The expressions (p3|), (p6|) depend only on the occupation probability Vi v ih = fihi^ — 



/s)(l — /k) of the 2p-2h states in the phase space. This fact leads to interpretation of the 
collision damping with linearized collision term as the relaxation process due to the coupling 
of one-particle and one-hole states to more complicated 2p — 2h configurations. 

The form of the collision integral fl2"S| ) in the Markovian limit {u — > 0) coincides with the 
standard expression for the collision integral in Fermi-liquid without retardation effects [|]. |3|] 
because in this case the term in square brackets of Eq.(EUl) tends to the value —S(Ae)/T. 



The equation (|23|) for some special explicit form of the quantity Xj was use d at first 



in Refs. |], ^T|, [22J. The derivation of the collision integral ( |2"5D is performed in Ref. JTB 
within framework of the extended time-dependent Hartree-Fock model. The expressions for 
the distortion functions Xj corresponding to a perturbation approach on collision term and 
including the amplitudes of the random phase approximation were used in this method. 
Note that the expression for the collision integral in two-component Fermi-system is 



obtained from Eq.(|26|) in the same manner as done in Ref.p3| under the assumption that 
chemical potentials and the equilibrium distribution functions are the same for protons and 
neutrons. 



3 Calculations of the relaxation times and nuclear mat- 
ter viscosity 

The collision integral can be used to calculate collisional relaxation times governing the 
dissipative properties of different physical quantities [@, |[ |§ |H| [□], [12], |26|. Be- 



low we calculate relaxation times, rj , of collective vibrations in two-component nuclear 
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matter consisting of neutron and protons subsystems. These collective relaxation times is 
determined by interparticle collisions within the distorted layers of the Fermi surface with 
multipolarity £: 

I />oo p /*oo r 

^ = J o de x J dn p J i t\p,e l )Y m {p)/ J q de, J dQ p 5f^Y m (p), (27) 

where Yi m (p) is the spherical harmonic function and J^'(p, e) are the linear combinations of 
the collision integrals for protons J p and neutrons J n in nuclear matter: J^> = (J p ± J n )/2; 
5f^ = (5f p ±5f n )/2. These times are proportional to the relaxation times defining the 
damping widths T^(L) of the isoscalar (+) and the isovector (— ) vibrations with frequency 
uj [12|, |23], |26], in regime of rare collisions with wr^' 3> 1 in the Fermi liquid. In particular, 
the collisional damping widths of giant resonances with dipole (L — 1) and quadrupole 
(L = 2) multipolarities resemble the widths in the relaxation rate approach 

T?\L)=%/t^\LI rt\L=l)=rtl ^(1 = 2)=^, (28) 

in the case when nuclear fluid dynamical model with relaxation is used |2B|, The colli- 
sional damping width [1^] of zero sound in the Fermi liquid with its relative velocity S r ~ 1 



is also given by Eq.(|28|) but with the use of the rj^^ oc T2 for t^(L). The time r^Zl 
at uj = is the thermal relaxation time determining the viscosity coefficient of the Fermi 
liquid [Q. 

The time t^I at uj = is the thermal relaxation time determining the viscosity coefficient 
of the Fermi liquid. 

The analytical expressions for collective relaxation times of the damping of the collective 
vibrations with frequency uj has the following general form in low-temperature and low- 



frequency limits (T, huj <C 6^)0, 123], E8 



-L = -L {{hujf + (2vrT) 2 } , -L = ^ [< a' av ^ +) > + < >] , (29) 

where a' av = (a' nn + <j' pp )/2; a'^, = dajj'/dfl is in- medium differential cross-section for scat- 
tering of the nucleons j and j' ( here, j = n or p, and similarly j' = p or n). The symbol 
< . . . > in Eq. (|2D|) denotes averaging over angles of the relative momenta of the colliding 
particles, 

<(...)>=- I* d(f)sm(4>/2) r d6{...). (30) 

7T JO JO 

The functions $^ define the angular constraint on nucleon scattering within the distorted 
layers of the Fermi surface with multipolarity £: 

$f } = 1 ± P,(cos0) - P e ((pspi)) T Pi((Pm)), (31) 
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where the scalar products {p3pi) and {pipi) are given as 

(p^i) = cos 2 (0/2) + sin 2 (0/2) cos 0, 
(p4Pi) = cos 2 (0/2) - sin 2 (0/2) cos 9. 



(32) 



Due to the momentum conservation and conditions pi= p F , the angle 9 agrees with the 
scattering angle in the center-of-mass reference frame of two nucleons. The angle defines 
the magnitudes of the relative momenta ki = {j>2 — Pi)/ 2 and kf = (p4 — pz)/2 before and 
after collision, respectively. The value of total momentum, P = pi + P2, also depends on a 
magnitude of the 0. We have 



kikf = k 2 cos 6, k 2 = k 2 = k)=p 2 F sin 2 (0/2), P 2 = Ap 2 F cos 2 (0/2). 



(33) 



Therefore the relative kinetic energy E re i of two nucleons as well as the energy of centrum 
mass motion E cm are dependent on angle 



E rel = k 2 /m = 2e F sin 2 (0/2), E cm = P 2 /2m = 2e F cos 2 (0/2) 



(34) 



and the total energy E tot = E re i + E cm holds only fixed, E tot = 2e F . Therefore the in-medium 
differential cross-sections a' i m of two nucleon scattering depend on the relative momenta ki 

and kf at fixed total energy rather then at fixed relative kinetic energy E re[ , because the 
magnitude of E re i changes with angle between colliding particles. The transfer momenta 
q = ki — kf = p 3 — pi and (f = —(ki + kf) = p\ — p 4 for scattering due to direct and 
exchange interactions respectively are also functions of the and 6: q = 2fc(0) sin(6 l /2) and 
q' = 2k((j>) cos(0/2). 

Now we estimate the collisional relaxation times in the case of the isotropic scattering 
with independent of energy the angle- integrated cross sections Ojy. Performing angular 



integration in (p9|) with the use of Eqs.(|30D and ( |3~1~D we find that 1/r 



(±) 



(±) 



and 



n 



„(±) 



(hu/27r) 2 + T 2 



a} 

2-(-lY 
2£ + l 



a 



m 

~2 



,(±) 

C£(Jav ~T CI £ CT n p 



(35) 



(-) 



1 -(-!)' 
21 + 1 



d 



(+) 



d 



(-) 



Q=o — Q=i 



0. 



where a av = [a pp +a nn +2a np ]/4 is the in-medium spin-isospin averaged nucleon- nucleon cross 
section. The magnitude of the in-medium cross section Ojy is taken usually proportional 
to the value of the free space cross section a^, with a factor F = &jj'/o'jf, so that the 
parameter can be rewritten in the form 



a 



(±) 



a^/F, af ] = 4.18/ c £ + 1.3d 



;(±) 



;(±) 



MeV. 



(36) 
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Here, the values a^J = 3.75 fm 2 and = 5 fm 2 are adopted || |27|; they correspond 
to the free space cross sections near Fermi energy. The dependence on I for relaxation times 
= j" 1 with the use of free space cross sections is shown in Fig. 1; huo = 13A3MeV . 
The Figs. 2,3 describe relative relaxation times versus multipolarity I for cross section 
parameterizations using Gogny and Skyrme effective two-body forces with parameters from 
[|15| , |30"1| . The relative relaxation times presented in Fig. 3 were calculated using cross section 



on Fermi surface. Solid and dashed lines connect the values which correspond to isoscalar 
and isovector modes of vibrations respectively. 

The magnitudes of the relaxation times are different for isoscalar and isovector modes of 
vibrations and they are dependent on the multipolarity I. The collisional relaxation times 
rather slowly vary with collective motion mode at isotropic scattering with energy inde- 
pendent free cross sections. As seen from the Fig. 3, the relaxation times calculated with 
effective interaction between nucleons are greater than that ones with cross-section in free 
space. It means that in-medium cross-section between nucleons in nuclear matter is smallest 
than in free space (~ 20 % for Skyrme forces and ~ 60 % for Gogny interaction). The relax- 
ation times Tg depend on frequency lu due to the memory effects in the collision integral. 
The temperature dependence arises from smearing out the equilibrium distribution function 
near the Fermi momentum in heated nuclei. The collisional rates 1/Vj are quadratic both 
in temperature and in frequency with the same relationship between the components much 
as in the zero sound attenuation factor of heated Fermi liquid within the Landau prescription 

©11,1 HI. 



According to response function approach the damping width of giant isovector dipole 



resonance (GDR) at the temperature T is presented in the form p9 



r ^ =2 " ( £?-5y + 1w (37) 

where 7 is determined by relaxation time r c {huo = E r ,T), 

7 = jz U w TV (38) 
T c (nuj = E r , 1 ) 

E r is an energy of giant dipole resonance; Eq = 41 A^ 1 ^ 3 MeV. Here, we determine quantity 
q from equality of the GDR width in cold nuclei with corresponding experimental value 
T exp : T(T = 0) = T exp The temperature dependence of the GDR width according to Eq.(|38|) 
in atomic nuclei 208 Pb and l20 Sn are shown on Figs. 6. Experimental data are taken from 



Refs. |pl|, |32| and they are indicated by points. The relaxation time r c is taken as equal 
to Tg = l: dash lines in Fig. 4 correspond to calculations with Skyrme forces; dash-dot lines 
- Gogny interaction; solid lines - calculations with free space cross-section for two-nucleon 
collisions. The temperature behaviour of the damping width is in rather close agreement 
with the ones of experimental data. 
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It is seen from Figs. 4 that in rare collision regime the dependence of the GDR widths on 
the collective vibration frequency and the temperature has the following form T oc (hu) 2 + 
4 (7rT) 2 , which corresponds to the Landau's prescription pD|, [531 [JJJ]. 

The relaxation times rather slowly vary with multipolarity of the Fermi surface dis- 
tortions governed by collective motion and two-body collisions. It gives possibility to use 
approximately the relaxation time ansatz for collision integral. 

This work is supported in part by the IAEA(Vienna) under contract 302-F4-UKR-11567. 
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Fig. 1: The relaxation times rff versus multipolarity I in cold nucleus 208 Pb with 
free space cross section. 
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Fig. 2: The relative relaxation times t^/t^ versus multipolarity I for different 
cross section parameterizations. 
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Fig. 3: The relative relaxation times Tg /r^f versus multipolarity i for different 
cross section parameterizations calculated on Fermi surface. 
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Fig. 4: The temperature dependence of the GDR width according to Eq.(|3~8|) in nuclei 208 Pb 
and 120 Sn: dash lines - calculations with Skyrme forces; dash-dot lines - Gogny interaction; 
solid lines - calculations with free space cross-section for two-nucleon collisions. 



